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Abstract 

We use the stochastic quantization method to obtain the free scalar propagator of a finite 
temperature field theory formulated in Minkowski spacetime. First we use the Markovian stochas- 
tic quantization approach to present the two-point function of the theory. Second, we assume a 
Langevin equation with a memory kernel and Einstein's relations with colored noise. The con- 
vergence of the stochastic processes in the asymptotic limit of the Markov parameter of these 
Markovian and non-Markovian Langevin equations for a free scalar theory is obtained. Our for- 
malism can be the starting point to discuss systems at finite temperature out of equilibrium. 
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1 Introduction 



The real time formalism allows one to discuss in detail finite temperature field theory out of 
equilibrium. The aim of this paper is to present an alternative approach to study finite temperature 
field theory in the the real time formalism using the stochastic quantization. 

The program of stochastic quantization, first proposed by Parisi and Wu [1] , and the stochastic 
regularization were carried out for systems described by fields defined in flat, Euclidean manifolds. 
A brief introduction to stochastic quantization can be found in Refs. [2] [3] [4], and a complete 
review is given in Ref. [5]. Recently Menezes and Svaiter [6] [7] implemented the stochastic 
quantization to study systems with complex valued path integral weights. Since we have a non- 
positive definite measure, the convergence of the stochastic process in the asymptotic limit of the 
Markov parameter is not achieved. To circumvent this problem, these authors assumed a Langevin 
equation with memory kernel and Einstein's relations with colored noise. In the asymptotic limit 
of the Markov parameter, the equilibrium solution of such Langevin equation was analyzed. It 
was shown that for a large class of elliptic non-Hermitean operators, which define different models 
in quantum field theory, the solution converges to the correct equilibrium state in the asymptotic 
limit of the Markov parameter r — > oo. 

We would like to remark that such kind of problems of non-convergence of the Langevin equa- 
tion in the stochastic quantization framework also appears if someone consider the stochastic 
quantization of classical fields defined in a generic curved manifold. For curved static manifolds, 
the implementation of the stochastic quantization is straightforward. In this situation it is possi- 
ble to perform a Wick rotation, i.e., analytically extend the pseudo-Riemannian manifold to the 
Riemannian domain without problem. Recently, the stochastic quantization for fields defined in 
a curved spacetime have been studied in the Refs. [8] [9]. Nevertheless, for non-static curved 
manifolds we have to extend the formalism beyond the Euclidean signature, i.e., to formulate 
the stochastic quantization in pseudo-Riemannian manifold, instead of formulating it in the Rie- 
mannian space, as was originally proposed. See for example the discussion presented by Hiiffel 
and Rumpf [10] and Gozzi [11]. In the first of these papers the authors proposed a modification 
of the original Parisi- Wu scheme, introducing a complex drift term in the Langevin equation, to 
implement the stochastic quantization in Minkowski spacetime. Gozzi studied the spectrum of 
the non-self-adjoint Fokker-Planck Hamiltonian to justify this program. See also the papers [12] 
[13]. Of course, these situations are special cases of ordinary Euclidean formulation for systems 
with complex actions. 

The main difference between the implementation of the stochastic quantization in Minkowski 
spacetime and in Euclidean space is the fact that in the latter case the approach to the equilibrium 
state is a stationary solution of the Fokker-Planck equation. In the Minkowski formulation, the 
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Hamiltonian is non-Hermitian and the eigenvalues of such Hamiltonian are in general complex. 
The real part of such eigenvalues are important to the asymptotic behavior at large Markov time, 
and the approach to the equilibrium is achieved only if we can show its semi-positive defmiteness. 
The crucial question is the following: what happens if the Langevin equation describes diffusion 
around a complex action? Some authors claim that it is possible to obtain meaningful results 
out of Langevin equation describing diffusion processes around a complex action. Parisi [14] and 
Klauder and Peterson [15] investigated the complex Langevin equation, where some numerical 
simulations in one-dimensional systems were presented. See also the papers [16] [17]. We would 
also like to mention the approach developed by Okamoto et al. [18] where the role of the kernel in 
the complex Langevin equation was studied. More recently, Guralnik and Pehlevan constructed 
an effective potential for the complex Langevin equation on a lattice [19]. These authors also 
investigated a complex Langevin equation and Dyson-Schwinger equations that appear in such 
situations [20]. 

We would like to remind that there are many examples in the literature where Euclidean 
action is complex. We have, for example, QCD with non-vanishing chemical potential at finite 
temperature; for SU(N) theories with N > 2, the fermion determinant becomes complex and 
also the effective action. Complex terms can also appear in the Langevin equation for fermions, 
but a suitable kernel can circumvent this problem [21] [22] [23]. Another important case that 
deserves attention is the stochastic quantization of topological field theories. The simplest case 
though is, of course, the stochastic quantization in Minkowski spacetime, as we discussed. This 
situation appears in the case for non-equilibrium problems, which are not amenable to an Euclidean 
formulation. Recently, Berges and Stamatescu [24] used stochastic quantization techniques to 
present lattice simulations of non-equilibrium quantum fields in Minkowskian spacetime. 

In the perturbation theory in quantum field theory at finite temperature there are three es- 
tablished methods: the Matsubara method [25] [26], the path ordered method [27] [28] and the 
Thermo Field Dynamics (TFD) approach [29] [30]. For non-equilibrium quantum field systems, 
the path ordered method or the Thermo Field Dynamics must be used. In the Thermo Field 
Dynamics approach one can develop finite temperature field theory in real time using the operator 
formalism, while in the path ordered method the path integral formalism is used. These three 
methods are related to each other by the analytical continuation of time variables. The motiva- 
tion of this paper is to present an alternative approach to study finite temperature field theory in 
the the real time formalism using the Markovian and the non-Markovian stochastic quantization 
procedures [31]. Basic ideas of the non-Markovian Langevin equation can be found in the Refs. 
[32] [33] [34] [35]. 

The outline of the paper is the following. Introduction is given in section I. In section II we 
present a brief review of the real time formalism in quantum field theory. In section III we use 
the Markovian stochastic quantization method to study a non-equilibrium thermal field theory 
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formulated in Minkowski spacetime. The non-Markovian approach of the stochastic quantization 
applied to a thermal scalar field theory is developed in section IV. Conclusions are given in section 
V. In the appendix, convergence conditions for the stochastic process are derived. In this paper 
we use h — c — ks — 1. 



2 Real-time formalism in finite temperature quantum field 
theory 

In this section we give a brief survey of the formulation of field theory at finite temperature in 
Minkowski spacetime. Unlike in the imaginary time formalism, in real time formulation, sums over 
Matsubara frequencies are absent and there is no need to analytically extend the Green functions 
back to Minkowski spacetime. Moreover, the real time formalism is the starting point for the 
development of the non-equilibrium quantum field theory, since the investigation of dynamical 
properties of systems is more naturally performed in this formalism. The real time formalism can 
describe non-equilibrium processes because the time variable plays a fundamental role and and 
cannot be traded in for an equilibrium temperature. 

For simplicity we work with a neutral scalar field. The field operator in the Heisenberg picture 
is given by 

0(t,x) = e^V(O,x)e-^, (1) 

where the time variable t is allowed to be complex. The main quantities to be computed are the 
thermal Green functions Gc(%i, ■ ■ ■ ,%n), defined as 

G c (x 1: ...,x N ) = (T c (0(xi) . . . <i>(x N ))) p , (2) 

where the time ordering is taken along a complex time path, yet to be defined. Considering a 
parametrization t = z(v) of the path, the following expressions: 

9 c (t-t') = 9(v-v% (3) 

5 c (t-t>) = \(v-v>), (4) 



define the generalized 9— and 5— functions. The functional differentiation is also extended in the 
following way: 

Ail t1 

5 c (t-lf)6*(x-x!), (5) 



Sj(x') 
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for functions j(x) defined on the path C. The Green functions defined by Eq. (2) can be obtained 
from a generating functional Zc[f3;j] through the expression 



Gc(xu ...,x N ) 



in the above equation, the generating functional is given by 



j=o 



Z[P;j] = Tr 



1311 T c exp J d 4 x j(x)<f>(x) 



(6) 



(7) 



= J V<P' {(j>'(x);t-iP\T c exp (i jf d 4 x j(x)<j>(xj) |0'(z); t), 

where the path C must go through all the arguments of the Green functions. It is also possible to 
note from this expression that the path C starts from a time ti — t and ends at a time tf — t — i/3. 
We may recast the generating functional into the form 



Z c [(3;j] = AT exp j-i J c d'x V }z£\/3;j] 

where Af is a normalization parameter and the free generating functional is given by 

Zg\P;j] = e X p{-±f c d i x j c d 4 yj(x)D^x-y)j(y)y 

In Eq. (9), the propagator Dq(x — y) is defined through the formula 

D£(x - x') = 9 c (t - t')D>(x, x') + 9 c (t' - t)D<(x, x'), 
where Dq(x,x') and Dq(x,x') are, respectively: 



D>(x,x') 
D<(x,x') 



(4>{x)(t>{x')) 
{4>{x')<t>{x))p. 



(8) 



(9) 



(10) 



(11) 



Since the propagator Dq(x — x') is properly defined in the interval —f3 < Im(t — t') < /3, one 
may conclude that the path considered must be such that the imaginary part of the time variable 
t is non-increasing when the parameter v increases. Furthermore, since we are interested in Green 
functions whose arguments are real, the path C must contain the real axis. One possible choice 
for the contour C is described in the following [36] : 



1. C starts from a real value large and negative. 
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2. The contour follows the real axis up to the large positive value — U. This part of C is denoted 
by d. 

3. The path from —U to —U — if, along a vertical straight line. This is denoted by C3. 

4. The path follows a horizontal line C 2 going from — 1 { — if to t { — if. 

5. Finally, the path follows a vertical line C4 from t, L — if to £j — i(3. 
Taking ti — > —00, the free generating functional can be factorized, 

Z^[(3;j] = Z^ iUC2 [(3;j]Z^ uC4 [(3;j\. (12) 

The Green functions with real time arguments can be deduced from Zq iUC2 [(3; j] only. The 
Zc :iUCi [P; j] generating functional can be considered a multiplicative constant. Choosing t and 
t' real, running from —00 to 00 and label the sources ji(x) = j(^ x ) an d h{ x ) — j(t — «/^/2, x). 

Also, one has ,° = 5 a b5 4 (x — x'). With this expressions one may rewrite the free generating 

0]b{X ) 

functional as 

Zg\P;j] = Af'exp |-± 1 d 4 :r | dV j a (z)L>^ 6 (*')}, (13) 

where, again, A/ 7 is a normalization parameter. The components of the matricial propagator 
D^ b (x — x') are given by 

D^x-x') = D F (t-t',x-x!), (14) 
D£ 2 (x-x') = D* F (t-t',x-x'), 
D[ 2 (x-x') = J D < (t-t' + i/3/2,x-x / ), 
D^(x-x') = D>(t-t' -i(3/2,x-x'). 

The effective generating functional can be written as 

Z c \J3-j] = J VfaVfc exp j-i J d 4 x d A x' <p a {x){D F l ) ab {x - x')<f> b (x')} (15) 



x 



exp j-i J d 4 x (V(<h) - V(<h)) +ijd 4 x j a (x)<j> a (x) 



The field 2 may be interpreted as a ghost field on the contour C 2 . This doubling of the field 
degrees of freedom, which does not occur in imaginary time formulation, is unavoidable in the real 
time formulation. For more details on this subject, the reader is referred to the original paper of 
Niemi and Semenoff [37] or the Landsmann and van Weert review 
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3 Real-time finite temperature quantum field theory: the 
Markovian stochastic quantization approach 

The real time formalism is a framework to describe both equilibrium and non-equilibrium systems. 
Dynamical questions, as for example a weakly interacting Bose gas having a temperature gradient 
can be studied only in the real time formalism, with a matrix structure of the propagator. Before we 
study the non-Markovian approach, in this section we will analyze the usual stochastic quantization 
of a finite temperature field theory formulated in Minkowski space. In Minkowski space, it is well 
known that the Langevin equation should be written as: 



— <f)(x,T)=l- 



+V(x,r), (16) 

4>{x)=4>{x,t) 



dr 6<f>(x) 
where S(4>) is the action for a free scalar field: 

S(<f>) = J d d x l - {dy - m 2 2 } , (17) 
and the correlation functions for the noise field are: 

(v(x,T))r, = 0, (18) 

( V (x,r) V (x',r')) v = 5(\r-r'\)5 d (x-x'). (19) 
If we consider a complex free scalar field, with an action S(4>, 4>*) given by 

S = J d d x {d^*d^ - m 2 <p*<p) , (20) 

we should have two Langevin equations, one for the scalar field and the other to its complex 

/ \ 

conjugate. If we write $ = , and working in Fourier space, we can write the Langevin 

\<f>*J 

equations as: 

^-$ a (k,T)=t(D^) ab <S> h (k,T)+r ]a (k,T), (21) 

{ V \ 

where we also consider a complex noise field, n — \ , a, b — 1, 2 and 

WJ 

, 1N , N f(k 2 -m 2 + ie) \ 
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The literature emphasizes [5] that the addition of a negative imaginary mass term — (i/2) e<p*(j) to 
the action in Eq. (20) is necessary in order to obtain convergence for the stochastic process being 
considered. So, it means that we can only take the limit e — > after all calculations have been 
performed. That explains the presence of the term — ie in the expression for the quantity (Dq 1 )^, 
that appears in Eq. (22). It is straightforward to obtain the two-point correlation functions 
for this case and to develop the perturbative solution to Eq. (21), with the following discussion 
on stochastic diagrams. We do not wish to go into details here. For the interested reader, we 
recommend the references [10] and [12]. 

Now, let us point our attentions to the non-equilibrium case. As we discussed in the last 
section, the doubling of the field degrees of freedom is unavoidable in the real time formulation. 

So, we can also write the field as an isovector = . The action for this isovector scalar field, 

in the free case, is given by: 

S = \j ' d A x A> a (i)( J D F 1 ) a() (x - x')(t> b {x'), (23) 

where the components of (D F ) ab are given by Eq. (14). In Fourier space: 

(D F ) ab (k) = (UXMiDoUimUUie), (24) 

where (D ) ab is the inverse of (Dq 1 )^, given by Eq. (22), and 

. . . . / cosh 9 sinh 6 \ , s 

M =Uh« (25) 

where: 

e P\ko\ 

cosh 2 # = -^-- -. (26) 

So, we can split (D F ) ab into two parts, (D F ) ab = (D ) ab + (Dp) ab , where (D ) ab (k) is temperature 
independent 

(DoU(k)=( 1 ^ JT l ° , ), (27) 
and all temperature dependence appears in (Dp) ab (k), which is given by 

-ie ( 2 sinh 2 9 sinh 26 \ 

(Dp) ab (k) = — _ m2)2 - ^ y s . nh2 ^ 2s . nh2e j. (28) 
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Notice that, in the limit e — > 0, we have: 

e 



(k 2 — m 2 ) 2 + e 2 



7r5(A; 2 -m 2 ). (29) 



As in the zero temperature case, we have to use the expressions above with finite e and take the 
e — > limit after all the calculations have been done in order to obtain convergence in the limit 
t — > oo. We also know that this is the case for the path integral formalism [37]. The Markovian 
Langevin equation for the non-equilibrium case is given by 

^Mk, r) = i(D^) ab (k)(f> b (k, r) + Va (k, r), (30) 

where the noise correlation functions are given by: 

( Va (k,r)) v = 0, (31) 
( Va (k,r) Vb (k',r')) v = (2n) d 5 ab 5 d (k + k')5(\r-r'\). (32) 

The solution for the Eq. (30) is given by: 

/oo 
dT'(g(k,T-T')) ab r] b (k,T'), (33) 

where g(k,r) = e lD p ^ t 9(t) is the Green function for the diffusion problem. In order to check 
convergence, i.e, to analyze if g(k,r)\ T ^, 00 — > 0, we must first diagonalize the matrix iDp l {k). 
From Eqs. (27) and (28), we have: 

2 sinh 2 9 „, o sinh 29 



where: 



Dp 1 (k) = I(k,e) [ k2 ~ m2 ~^ (fc 2 -™ 2 ) 2 +^ 2 (fe 2 -m 2 ) 2 +e 2 I (34) 

\ (t ! -m ! )V ™ " fc 2 -m 2 +i e (fc2_ m 2 )2+e 2 Z MM (7/ 

-(A; 2 - m 2 ) 2 - e 2 (cosh 4 + sinh 4 0) - f sinh 2 26 

I{K e) ~ ((jfc2- m 2)2 + £ 2) • ( 35 ) 



Diagonalizing iD F 1 (k), we get the matrix D'(k), given by: 



^(fc) =</(*,£) ( A J ° ), (36) 



where: 



±\]{k 2 - m 2 ) 2 - e 2 sinh 2 2# - % e(l + 2 sinh 2 0) 
A ± = (A; 2 - m 2 ) 2 + e 2 ' (3?) 



Since I(k,e) < 0, we notice from the above equations that, indeed, we get g(k,r)\ T ^ 00 — > 0. We 
also remark that, as in the zero temperature case, the convergence of the stochastic process was 
possible because we have maintained in the Eqs. (27) and (28) a finite e. As the reader can easily 
verify from the Eqs. (36) and (37), if we take the e — > limit in the beginning of the calculations, 
we should lose the convergence factor e - / ( fc . £ M 1 + 2sinh e ), g Q this limit should be taken after all the 
calculations have been done in order to obtain convergence in the limit r — > oo. We also know 
that this is the case for the path integral formalism [37]. 

Now, we are ready to calculate the two point function ((f) a (k, r)4>b(k', t)) v . Proceeding with sim- 
ilar calculations as the zero temperature case, it is possible to show that the two-point correlation 
function is given by: 

(Uk, r)Mk', r)) v = (2n) d 5 d (k + k')i (D F ) ac (k)(l - e 2iD ^^ T ) cb , (38) 

so we see that, in the limit r — > oo, we recover the usual result. We are interested now to see the 
effects of a memory kernel in this non-equilibrium quantum field theory. This is the subject of the 
next section. 



4 Real-time finite temperature quantum field theory: the 
non-Markovian stochastic quantization approach 

The aim of this section is to study finite temperature quantum field theory in Minkowski spacetime, 
using the non-Markovian stochastic quantization approach. In Minkowski space, the Langevin 
equation with memory kernel is written as 



d , , , . r . . , , , ss 



—<j)(x,T)=i r ds m a ( 

OT JO 



5(j)(x) 



+ V(x,r), (39) 

4>(x)=<t>{x,s) 



where S is the action for the free scalar field, given by Eq. (23). The noise field distribution is 
such that its first and second momenta are given by 

( V a(x,r)) v = 0, (40) 
( Va (x,r) Vb (x',r')) v = 25 ab M A (\r-r'\)5 d (x-x'), (41) 

that is, the distribution is a colored noise Gaussian distribution. We remind the reader that Eq. 
(39) is to be understood as a matrix equation. 
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Using a Fourier decomposition for the scalar and noise fields, given by 

X(x,t) = -±-j [d d k e ikx X(k,r), 



(42) 



where the field X represents either the noise field r] and the scalar field 0, we obtain that each 
Fourier mode (f>(k, r) satisfies a Langevin equation of the form 



d_ 



(f) a (k, t) =i [ ds M A (r - s) (D F l ) ab (k) (f) b (k, s) + rj a (k, r), 

JO 



(43) 



where D F 1 (k) is the inverse of Dp(k), defined by Eq. (24). With this decomposition, we obtain 
from Eq. (40) and Eq. (41) the following relations for the noise field Fourier components, 

(Va(k,r)) v = 0, (44) 
(r j a(k,T)r ib (k>,T')) r , = 2(2n) d 6 ab M A (\r-r'\)5 d (k + k'). (45) 

Defining the Laplace transform of the memory kernel as 

M(z) = / dr M A (r)e- 2T , (46) 
Jo 

we obtain the solution for Eq. (43), subject to the initial condition 4> a (k,r = 0) = 0, a = 1,2: 

roo 

<f> a (k, r) = / dr' G ab (k, r - r')r] a (k, r') (47) 
Jo 

In Eq. (47) G ab (k, r — t') = fl ab (k, r — t')0(t — t') is the retarded Green function for the diffusion 
problem and each component of the fi-matrix is defined through its Laplace transform, 

-l 



Qab — 



z5 ab -iM{z)(D F l ) (k) 

ab 



(48) 



Thus, the two point correlation function in the Fourier representation is written as 
(Mk,r)Mk',r')) ri = D ab (k,r,r') 



= 2(2tt) d 5 d (k + k') J o ds J o ds' [Q(k,r - r')Q(k,r - s')] ab M A (\s - s'\). 
The two-dimensional Laplace transform of the above equation is given by 

roo roc pt rr' 

/ dr e- ZT / dr' e- ZT / ds / ds' [Q(k, r - r')tt(k, r - s')] nb M A (\s - s'\ 
Jo Jo Jo Jo 

' 'M(z) + M{z')\ 



(49) 



[Q{k,z)Q{k',z')} 



ab 



z + z' 



(50) 



10 



Using the Eq. (48) this expression becomes 

roc roc rr rr' 

/ dre- ZT / dr' e~ ZT / ds / ds' [Q(k, r — r')Q(k, r — s')] ah M A (\s — s'\) 
Jo Jo Jo Jo 

. fn(k,z) + n(k,z') 



= I 



n(k,z)n(k,z')\ (D F ) cb (k). (5i) 



z + z< 

Applying the inverse transform, we obtain for the two-point function 

D ab (k, r, t') = 2i(2v) d 5 d (k + k') (Q(k, \r' - r|) - Q(k, r)Q(k, r')) ac {D F ) ch {k). (52) 

In order to investigate the convergence of the above equation, we need to specify an expression 
for the memory kernel Ma- We set 

M A (r) = ±A 2 e - A X (53) 

Substituting the Laplace transform of the Eq. (53) in the Eq. (48), we have that the fi- matrix 
is given by 

/n u (fc,r) n 12 (k,r)\ 
to(k,T) = , (54) 

\il 21 (k,T) U 22 {k,T)J 

where the components VL ab (k, r) are given in the Appendix. So, we are in a position to present an 
expression for the two-point correlation function in the limit r = r' — > oo: 

D a b(k, r, r^lr^/^oo = ^ (27r) d 5 d (&; + k')(D F ) ab (k), (55) 

so that, in the limit e — >■ 0, we have: 

£>ab(A;, r,r') Ir^ooj^o = i (27r) d <5 d (£; + k')(D F ) ab {k)\ e ^ - (56) 

A question still remains opened. What are, if any, the advantages of our non-Markovian method 
over the usual Markovian one? In order to answer such question, we shall apply a Fokker-Planck 
analysis. As we know, correlation functions are introduced as averages over rj: 

( <f>(xi, r 1 )<f>(x 2 , r 2 ) ■ ■ ■ <f)(x n , T n ) )„ = 

N J [drj\exp(--J d d x J dr r] 2 (x, r)) ^>(x u n)(f)(x 2 , r 2 ) • • • (f)(x n , r n ), (57) 
where obeys Eq. (16) and N is given by: 

AT- 1 = J [d v ] exp(-^ J d d x J drr] 2 (x,r)^ (58) 
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An alternative way to write this average is to introduce the probability density P[cj>, r], which 
is defined as [39]: 

P[0,r] = J [dT,]exp(-±J 'd d x J dr V 2 (x,r^ J] Wv) - (f>(y, r)). (59) 

In terms of P, the correlation functions will read: 

( (f)(x 1 , n)(f)(x2, t 2 ) • • • 0(x n , r„) ) v = N [d(f)} 4>(xi, T 1 )(f)(x2, r 2 ) • • • 0(x„, r n )P[0, r]. (60) 



The free probability density P satisfies the following Fokker-Planck equation: 

^ I dd *ih (ik ~ m T] • (61) 

where S is given by Eq. (17) and with the initial condition: 

^0] = IIW- (62) 
y 

The stochastic quantization says that we shall have: 

i- du i exp(i,S[0]) 
w hm P 0,r = — — 63 

/[d0]exp(«S'[0]) 

where the limit is supposed to be taken "weakly" in the sense of the reference [39]. 

In our real time non-Markovian case, if we notice the resemblance between our retarded Green 
function G a b(k, r) and the one found in Ref. [31], we may follow similar steps to calculate the free 
probability density. It is given by, in momentum space: 

P[0,r] = TV^exp Q J dk (j) a (k)D~^(k, r, r)(f) b (-k)j (64) 

where iV _1 is a normalization factor and D~^(k, r, r') is the inverse of D ab (k, r, r'), defined by Eq. 
(52). It is easy to verify that, in the limit r — > oo, P[</>,r] will satisfy, up to constants, similar 
relations as obtained by Ref. [39]. However, for massless scalar theories, those estimations in 
such reference decay as inverse power of r. In our approach, even in the massless situation, an 
exponential behavior is found. Therefore, in the limit r — > oo, it seems that we get an improved 
convergence. 

In conclusion, we have used the stochastic quantization method to study thermal field theory 
formulated in Minkowski spacetime. We have assumed a Langevin equation with a memory kernel 
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and Einstein's relations with colored noise. From the above last equation, we see that the equilib- 
rium solution in the asymptotic Markov time of this non-Markovian Langevin equation for scalar 
theories can be obtained. Our approach based in stochastic quantization using a non-Markovian 
Langevin equation proved to be well suited to quantize a classical field out of equilibrium in the 
real time formalism at finite temperature. 



5 Conclusions 



In the past several years there have been a lot of interest in quantum field theory at finite 
temperature. There are three main formulations of finite temperature field theory. The imaginary 
time approach or the Matsubara formalism and the real time formalism, which can be operatorial 
or use the path integral approach. In the real time formalism, it is necessary to double the number 
of the field degrees of freedom. To quantize a classical thermal field theory out of equilibrium, using 
the stochastic quantization, we are forced to work in the Minkowski spacetime, where naturally 
a imaginary drift term appears in the Langevin equation. Since in this case the path integral 
weight is not positive definite, the stochastic quantization in this situation is problematic. Parisi 
and Klauder proposed complex Langevin equations [14] [15] , and some problems of this approach 
are the following. First of all, complex Langevin simulations do not converge to a stationary 
distribution in many situations. Besides, if it does, it may converge to many different stationary 
distributions. The complex Langevin equation also appears when the original method proposed by 
Parisi and Wu is extended to include theories with fermions [21] [22] [23]. The first question that 
appears in this context is if make sense the Brownian problem with anticommutating numbers. 
It can be shown that, for massless fermionic fields, there will not be a convergence factor after 
integrating the Markovian Langevin equation. Therefore the equilibrium is not reached. One way 
of avoiding this problem is to introduce a kernel in the Langevin equation describing the evolution 
of two Grassmannian fields. 

In this paper, we have used the method of the stochastic quantization to study thermal field 
theory formulated in real time. As we discussed, this closed time path method can be used to 
describe non-equilibrium thermal field theory. First we use the Markovian stochastic quantization 
approach to present the two-point function of the theory. Second, we assumed a Langevin equation 
with a memory kernel and Einstein's relation with colored noise. The equilibrium solution of such 
Langevin equation was analyzed. We have shown that for a large class of elliptic non-Hermitean 
operators which define different models in quantum field theory converges in the asymptotic limit 
of the Markov parameter r — > oo, and we have obtained the free Green functions of the theory. Al- 
though non-trivial, the method proposed can be extended to interacting field theory with complex 
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actions, where a consistent perturbation theory out of equilibrium can be developed. 

Finally, the literature has emphasized that the stochastic quantization is only an alternative 
formalism to quantize a classical field theory, but new results have not been obtained. Neverthe- 
less the stochastic quantization and the Langevin equation can be extremely useful in numerical 
simulations of field theory models [40] [41]. The implementation of this non-Markovian Langevin 
equation on the lattice is under investigation by the authors. 
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A Appendix 

In this appendix, we derive the components fl a b(k, r). We can write the Q- matrix as an inverse of 
the matrix A, whose components are given by 



A _ fz-iM(z)d' iM(z)b' 
' ~ V iM{z)b' z-iM(z)a' 

The quantities that appear in the A-matrix are defined by 



and 



(A.l) 



1 ie 2 

k 2 — m 2 + ie {k 2 — m 2 ) 2 + e 2 ^ ^ 

— if 

sinh20, (A.6) 



[k 2 — m 2 ) 2 + e 2 
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So, we will have 

{Q) ab {k,z) = [A )ab=[ n f , v n f , J, (A.8) 

\tt 2 i{k,z) il 2 2{k,z)J 

where: 

, . 2;-iM(z)a' 
ll( ,Z} = (z-iM(z)a')(z-i M{z) d') + M 2 (z)b> 2 ' ( } 

n^MHM*) = ^ _ . M ^ - /} ^ '_^ ( * } d /} - M2 ^ , (A.io) 

, , z-iM(z)d' 

M ,Z) ~ (z-iM(z)a')(z-iM(z)d') + M*(z)b"' { } 

The Laplace transform for the memory kernel, Eq. (53), is given by 

"W = TITiF < a ' 12 > 

So, inserting this result in Eqs. (A. 9), (A.IO) and (A. 11), we get: 

MM = (A.13) 

n 12 (k,z) = n 21 (z) = -^ Z + ^ A2 \ (A . 14) 

where tj = i 3 — 2 — , j = a, 6, d, and 

= z 3 + 2A 2 z 2 + (A 4 -t^z-tj A 2 , (A.16) 
Q(s) = z 4 + 2A 2 z 3 + (A 4 -u)z 2 -uA 2 z + v, (A.17) 

with u = ii^+pAl an d ^ = ^ffl^ , From Eqs. (A.2), (A.3) and (A.4), we have that: 

_ A 2 e(l + 2 sinh 2 6)((k 2 - m 2 ) 2 + e 2 ) 

-(P -m 2 ) 2 -e 2 (cosh 4 + sinh 4/ 9) + f sinh 2 20' 
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and 



A 4 



(( jfe 2_ m 2 ) 2 + e 2 ) 



4 -(P-m 2 ) 2 -e 2 (cosh 4 # + sinh 4 #) + f sinh 2 2#' 

so u < and t> > 0. In order to get the inverse Laplace transform of each component of the 
f2-matrix, we must seek for the solutions of the quartic equation Q(z) = 0. As it is well known, a 
general quartic equation is a fourth-order polynomial equation of the form 



z 4 + (I3Z 3 + a 2 z 2 + a\z + ao = 0. 



(A.20) 



Using the familiar algebraic technique [42], it is easy to show that the roots of Eq.(A.20) are given 
by: 

--03 + -R+-D, 



where: 



Z\ 
Z2 
Z 3 = 
Za = 



1 In 1 ^ 

- A a 3 + -R--D, 
1 1 1 



\r-\e, 

2 2 



7(13 



n f 1 2 V /2 



D = < 



x 1/2 



(F(R) + Gj for R ^ 



1/2 



F(0) + H) for R = 0, 



E = I 



F(R)-G) forR^O 

x 1/2 

F(0) - for i? = 0, 



(A.21) 
(A.22) 
(A.23) 
(A.24) 

(A.25) 
(A.26) 



(A.27) 



F(R) = -a 2 3 



H = 2[yi-4a 



R 2 - 2a 2 , 
1/2 



(A.28) 
(A.29) 
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G= i(4a 3 a 2 -8ai - dj)ir\ (A.30) 
and yi is a real root of the following cubic equation: 

V 3 - d2y 2 + (aia 3 - 4a )y + (4a 2 a -a\ — a 2 a ) = 0. (A.31) 

For convenience, let us assume that R, defined by Eq.(A.25), does not vanish. Comparing Eqs. 
(A. 17) and (A. 20), we easily see that a 3 = 2 A 2 , a 2 = A 4 — u, a\ = -uA 2 and do = v. Therefore, 
the inverse Laplace transform of Q ab is given by: 

nn(*,r) = S {zi,ta l s eZ1T + 

(z 1 - z 2 )(z 1 - z 3 ){z l - z A ) 

P(z 2 ,t a ) ^ Z2T 



e" 2 ' + 



{z 2 ~ z 1 )(z 2 - z 3 )(z 2 - z A ) 



P(z3,t a ) 



(z 3 - zi)(z 3 - z 2 )(z 3 - Z A ) 
P(z A ,t a ) 

(Z4 ~ Z^){z A - Z 2 )(z A - Z 3 ) 



e Z4T , (A.32) 



n 12 (fc,r) = n 21 (*,r) = -L tb * 1 + tb u S ^ T + 

\{zi - z 2 ){zi - z 3 )(zi - z A ) 
I t b z 2 + t b A* + 

(Z 2 ~ Zi)(z 2 - Z 3 )(Z 2 ~ Z A ) 

, t b z 3 + t b A 2 ^ + 

(z 3 - zi)(z 3 - z 2 )(z 3 - Z A ) 

^ii±M? >A (A .33) 

{z A - Zi){Z A - z 2 )(z A - z 3 ) j 

and, finally, fl 22 (k, r) = Qn(k, r; t a — > to)- The roots Zi are given by: 

A 2 1 1 , kni . 
z 1 = -— + -ia + -i lf , (A.34) 

A 2 1 1 . . 

z 2 = -— + -ta - -i 7 , (A.35) 

A 2 1 1 . . 

z 3 = -— - -ta + -i 7 , (A.36) 
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A 2 1 1 

z± = -— - -ia - -vy, (A.37) 

with a = (\u\ — yi) 1 / 2 and 7 = (—A 4 + \u\ + yi) 1 ^ 2 being real quantities. The Eqs. (A. 32) and 
(A. 33) can be rewritten in a simpler form as: 

n n(k , T) _ _( H (^ T ) + _^_ sin (0L±l) T )) ftl + 



+8t „sin(^) S in(^) + 



+4«„A 2 ( JlS m(^±^r) + 92S in(^— ^t) ; I-— , (A.38) 



-2 sin 




where 



/ii = -(a+ 7 ) 2 -A 4 , (A.40) 

^2 = (a- 7 ) 2 + A 4 , (A.41) 
2A 2 

91 = i-i-rr^' (A - 42) 

(cr+ 7) 
2A 2 

32 = i- 7 r, (A.43) 

(a - 7) 

and, as before, f222(^ ? T) = T;t a — > t^). Let us consider the convergence of the stochastic 

process. In order for our stochastic process to converge, the retarded Green function for the 
diffusion problem should obey G a b(k, t)\ t ^oo — > 0. In other words, we must have fi a fe(^, t)|t^oo — ► 
0. From these last expressions, it is easy to see that the stochastic process will converge, if the 
quantities a and 7 are real, as imposed before. This lead us to the following conditions: \u\ —y\ > 
and |u| + y\ — A 4 > 0, or, combining those requirements, \u\ > 4r- Remembering Eq. (A. 18), we 
will have the following convergence criterium: 

e(l + 2sinh 2 #)((A; 2 -m 2 ) 2 + e 2 ) A 2 

> — , (A.44) 



(k 2 - m 2 ) 2 + e 2 (cosh 4 9 + sinh 4 9) - f sinh 2 29 2 ' 
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Now let us present the quantity y\. As was stated before, y\ is a real root of a cubic equation: 

z 3 + b 2 z 2 + b 1 z + b = 0. (A.45) 

Comparing Eqs. (A. 17), (A. 20), (A. 31) and (A.45), we have the following identifications: b 2 = 
u - A 4 , bi = -2 A 4 u - 4 v and b = -4 u v - w 2 A 4 . If we let: 

l b 1 -hl (A.46) 



3 9 

-(6162- - 
6 V 7 27 



\b 1 b 2 -3b )-±-b 3 , (A.47) 



we will have that 



4 A 4 A 8 4 w 2 /a 

9 = -AS, — -5.-5, (A.48) 

r = iAS u + lA* v + ±A* u * + luv-±A" + ±u*. (A.49) 

So, writing si = (r + \/g 3 + r 2 ) 1 / 2 and s 2 = — Vl 3 + J" 2 ) 1 ^ 2 , we have that: 

A 4 - u 

Vi = (si + s 2 ) + — - — . (A.50) 
As one can see, y\ > 0. Also, in the limit e — > 0, y\ becomes a polynomial of A. 
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